Can MONDian vector theories explain the cosmic speed up ? 
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Generalized Einstein - Aether vector field models have been shown to provide, in the weak field 
regime, modifications to gravity which can be reconciled with the successfull MOND proposal. Very 
little is known, however, on the function J-(JC) defining the vector field Lagrangian so that an 
analysis of the viability of such theories at the cosmological scales has never been performed. As 
a first step along this route, we rely on the relation between J-(K.) and the MOND interpolating 
function fi(a/ao) to assign the vector field Lagrangian thus obtaining what we refer to as MONDian 
vector models. Since they are able by construction to recover the MOND successes on galaxy scales, 
we investigate whether they can also drive the observed accelerated expansion by fitting the models 
to the Type la Supernovae data. Should be this the case, we have a unified framework where both 
dark energy and dark matter can be seen as different manifestations of a single vector field, ft 
turns out that both MONDian vector models are able to well fit the low redshift data on Type la 
Supernovae, while some tension could be present in the high z regime. 

PACS numbers: 98.80.-k, 98.80.Es, 95.36.+d, 95.36.+X 
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I. INTRODUCTION 

The old standard cosmological model of a universe full 
of baryon matter only and regulated by the laws of Ein- 
stein General Relativity provided a remarkable picture to 
consistently explain the background evolution of the uni- 
verse from the initial singularity up to the present day. 
Unfortunately, two observational facts then putted in se- 
rious trouble this elegant scenario. On the one hand, 
it turned out that the gravitational field measured on 
different scales may not be reconciled with what is in- 
ferred from baryonic matter only, the most famous ex- 
ample being represented by the flat rotation curves of 
spiral galaxies On the other hand, recent data on the 
Hubble diagram of Type la Supernovae (hereafter SNela) 
provided the first evidence for an accelerated expansion 
[2|] . Surprising as it was, the cosmic speed up has been 
then confirmed by updated SNela data 0, EL IE H 0, 
the anisotropy spectrum of the Cosmic Microwave Back- 
ground Radiation (CMBR) as measured by both balloon 
Q and satellites 0, [l(| experiments, and the data on the 
large scale clustering of galaxies * J. jj as estimated from 
large spectroscopic galaxy surveys. 

Both these problems have been traditionally solved by 
adding new ingredients to the cosmic pie, while leav- 
ing unchanged the underlying theory. Cold dark matter 
(CDM) has been then invoked to fill the gap between the 
baryon mass (constrained by the Big Bang nucleosynthe- 
sis) and what is needed to reproduce the observed gravi- 
tational field. Moreover, this mysterious component must 
not interact with anything but gravitationally so to have 
escaped any detection notwithstanding the incredible ef- 
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forts spent up to now to find direct evidence for any dark 
particle. At the cosmological scales, dark matter still be- 
haves as normal matter so that a further component is 
needed in order to drive the accelerated expansion. This 
new actor on the scene must, moreover, have an unusual 
negative pressure and not cluster on galactic and clus- 
ter scales in order to not oppose gravitational clustering. 
Referred to as dark energy, the nature and nurture of this 
new term is still debated with the cosmological constant 
A [l2l] being the simplest candidate. Although added by 
hand ad hoc, the cosmological model comprising both 
of them and known as ACDM is able to fit extremely 
well the full dataset at disposal [13| thus worthing the 
name of concordace model. Notwithstanding this remark- 
able success, the ACDM is theoretically unappealing be- 
cause of several well known shortcomings which moti- 
vated the search for other dark energy candidates such 
as quintessence scalar fields 14| and modifiedgravity the- 
ories, either introducing higher dimensions [15j | or corret- 
ing the gravitational Lagrangian 

nana. 

It is worth noting that both dark matter and dark 
energy can be avoided if one accepts that the problem 
is not with what is missing, but rather with how we 
describe gravity. From this point of view, all the evi- 
dences for dark matter and dark energy should be rather 
seen as evidence that something is wrong with the New- 
ton-Einstein theory of gravity. As a next step, one has 
therefore to look for a way to modify gravity at both the 
galactic and cosmological scales possibly finding a unified 
mechanism explaining phenomena taking place on very 
different scales such as the flatness of rotation curves and 
the acceleration of the universe expansion. 

Generalized Einstein- Aether theories are one of these 
proposed mechanisms based on the introduction of dy- 
namical timelike vector field with noncanonical kinetic 
terms. Such a proposal is not fully new, but builds upon 
the extensive analysis of standard Einstein - Aether theo- 
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ries |l8[ as phenomenological probes of Lorenz violation 
in quantum gravity. Indeed, the nonvanishing expec- 
tation value of the Aether field will dynamically select 
a preferred frame in the spacetime, namely the one in 
which the time coordinate basis vector dt is aligned with 
the direction of the vector field, thus leading to viola- 
tion of the Lorenz and gauge invariance. The interest in 
these models was then renewed after the demonstration 
[20T | that TeVeS, i.e. the fully relativistic and complete 
MOND theory , can be reformulated as a generalized 
Einstein - Aether theory. Later, it has been shown that 
it is possible to get a MOND - like behaviour in the weak 
field regime of vector models [2lj. Motivated by these 
considerations, we explore here the cosmological viabil- 
ity of these theories by suitably setting the vector field 
Lagrangian in such a way that the successfull MOND 
phenomenology is recovered in the Newtonian limit. 

The structure of the paper is as follows. The basic 
equations and properties of the vector models arc de- 
scribed in Sect. II where we also explain how we choose 
the vector field Lagragian. Sect. Ill is devoted to a de- 
tailed description of the statistical methodology and the 
data used to test the models in the low redhift regime 
with a discussion of the results obtained from the like- 
lihood analysis. A similar test is then presented and 
discussed in Sect. IV where we add high redhift data to 
probe this complementary redshift range. Conclusions 
and perspectives are finally given in Sect. V. 



II. THE MODEL 

A general action for a vector field A coupled to gravity 
can be written as 
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where Gn is the Newton gravitational constant, g the 
metric (with the signature -, +, +, +), R the Ricci scalar 
and Sm the matter action. The vector field Lagrangian 
C may be a whatever covariant and local function, but, 
following [2l[ , we will only consider the case : 
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Lagrangian through the generic function T(JC) thus al- 
lowing to have also a noncanonical contribution by these 
terms. The Ci quantities are dimensionless constants, 
while M is a scaling mass parameter and A a nondynami- 
cal Lagrange multiplier with dimensions of mass squared. 

It is possible to show [HI, [13, HI that the Einstein field 
equations may still be formally written as 



G a p — T a p + 8ttGnT^ 



(5) 



with G a p and the usual Einstein and matter stress - 

energy tensors, while T a p contains the terms related to 
the vector field and some of its derivatives (see, e.g., [2l[ 
for its full expression), while the equation of motion for 
the vector field reads : 



with T' — dT /dJC and 
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and defined the functional derivative : 
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Introducing the flat Robertson - Walker metric in Eq. (J5]) 
then gives the modified Friedmann equations which read : 



1_oKV2 « 



d 
dt 
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= 8t:G n (p + p), (10) 
where, because of the metric symmetry, it is : 



JC = 3a 



H 2 3a H 2 



M 2 



e 2 m 



(11) 



Here, we have defined two new parameters, namely e — 
M/Hq and a = c\ + 3c 2 + C3, while (p,p) are the total 
energy density and pressure of the source terms (matter, 
radiation, neutrinos, etc.) and H = a/a is the usual 
Hubble parameter. Hereafter, we will denote derivatives 
with respect to t with a dot and to K, with a prime, and 
use a subscript "0" to label present day quantities. 

It is convenient to rearrange the above equations in a 
different form. To this end, we first solve Eq.Q to get : 



It is worth noticing that K, represents all the possible 
canonical kinetic terms that we introduce in the vector 
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with E(z) — H(z)/Hq the dimensionless Hubble param- 
eter. Inserting Eq. (fl"!?|) into Eq. pu]) and using 



d 
dt 



(1 + z)H E(z) 



dz 



to change variable, one finally gets : 



p(a/a ) = /i(v K) = 2 + ci 



dT 
dIC 



(18) 



Note that, by this position, the MOND acceleration scale 
ao turns out to be related to the mass scale M as 



d 

dz 



e 2 T(K) P / Pc 



6E 2 (z) E 2 (z) 
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E\z) = 



2>{p + p) / Pcrit 



(l + z)E(z) 

_ (13) 

with p cr it — 3Hq /8itGn the present day critical density. 
Let us now assume that the universe is filled by a matter - 
like term and radiation. Since the vector field is only 
coupled to gravity, the continuity equation 

p + 3H(p + p) = 

still holds for both matter and radiation and we can 
therefore write the rhs of Eq. (fL3|) as 



(14) 



Z(P + P)/ Peru = [3» M + 4q.(l + z)](l + z) 2 
(l + z)E(z) 



E(z) 



with f2j = pi(z — 0)/ p cr n the present day density pa- 
rameter of the i-th component. 

Eqs. (fl2|) - (p~3|) clearly show that a key role in deter- 
mining the cosmic evolution is played by the functional 
expression adopted for T{JC). A guide to the choice of 
this quantity is provided by the observation that, in the 
nonrelativistic regime, the field equations reduce to : 
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with 



K, = —c\ 
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being $ the gravitational potential sourced by the density 
distribution p. By using the above expression, we get : 
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where we have reintroduced the speed of light c. It is 
worth noting that such a model thus provides a natural 
mechanism to explain why one observationally finds that 
ao ~ cHq which now emerges as a consequence of the 
local and cosmological phenomena being different mani- 
festations of the same underlying theory. 

Since we know that MOND makes it possible to fit the 
flat rotation curves of spiral galaxies (see, e.g., [24[ and 
refs. therein), it is reasonable to assume that a viable 
expression for T(K,) should lead to the same interpolat- 
ing function p(a/ao) which is successfully used on local 
scales. Two such functions are the simple form [25| 



p(x) 



and the standard form [2t 
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Although Eq. (l2"Tj) has been the first proposal to be tested 
with success [24j , recent analyses [25|, [27| seem to favour 
Eq. (|20|) . However, in order to gain further insight on 
the problem of which interpolating function is better 
motivated, we will consider both cases contrasting them 
against data probing radically different scales. 

Inserting alternatively Eas.(|20]) and (J^TJ) into Eq.([Ig]) 
gives respectively : 



T{jq = — 

ci 



n/£ - In ( 1 + V£ 



(22) 



V • M|V$|/M)V$] = 8wG N p 



(17) 



which is the modified Poisson equation for the MOND 
theory 1 , provided we identify the MOND interpolating 
function p(a/ao) with 



1 Actually, such a modified Poisson equation was proposed for 
MOND in the framework of the so called AQUAL theory [2^1 . 
one of the first attempts to work out a relativistic MOND theory. 
AQUAL was later abandoned since it turned out to be unable to 
account for lensing data without cold dark matter thus being in 
contrast with the original MOND philosophy. 



F(JC) = \K- y/K{l + K) + In (VJC + VT+K, 



2 

ci 

(23) 

Hereafter, we will refer to the models assigned by 
Eqs. ([22|) and ([23]) as the simple and standard MONDian 
vector model in order to clearly remember that the main 
ingredient of this kind of cosmological theory is the pres- 
ence of a vector field with a Lagrangian constructed in 
such a way to recover the simple or standard MOND in- 
terpolating function. Both expressions for !F(K) depend 
on \/lC so that it must be a > in the cosmological 
setting. In the small K. limit, T{K) ~ K, for both the 
simple and standard MONDian vector models so that we 
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can use the perturbative analysis in [28| for the Einstein - 
Aether models to see whether the condition a > is com- 
patible with the constraints on the coefficients c,. The 
requirement that the Hamiltonian for the perturbations 
is positive definite implies c\ < 0, while the constraint 
(ci + C2 + ca)/ci > has to be set in order to avoid tachy- 
onic propagation of spin - modes. On the other hand, if 
we allow superluminal propagations of both spin - and 
spin -2 modes, as supported in [29(, we then get that for 
c 2 > and 

-(c! + 3c 2 ) < c 3 < -(d + c 2 ) 



which can be solved numerically for given values of 
(f2jvf> fir> h). Actually, since Q r is typically set by the 
CMB temperature, the simple MONDian vector model 
is fully characterized by only two parameters. This is 
the same as the concordance ACDM model where the 
same parameters (QM,h) have to be assigned in order 
to compare the theory with the data. However, in that 
scenario, two different ingredients are invoked in order to 
explain the dynamics of galaxies and the cosmic speed 
up, while here the solution to both problems comes out 
as a consequence of the presence of a single vector field. 



it is indeed a > so that we can safely consider 
the MONDian vector models without violating any con- 
straint on the Ci coefficients. 



A. The simple MONDian vector model 



B. The standard MONDian vector model 

We can repeat the same steps as before to get the mas- 
ter equation for the case when T{1C) is given by Eq. (f2"3")l . 
Not surprisingly given the similarities of the models, we 
still get an algebraic relation : 



In order to determine the dynamics for this case, we 
have to insert Eg. ([22)1 into the master equation (JT3J) and 
solve with respect to E(z). Somewhat surprisingly, after 
some algebra, we get : 

~ r ~ / , n dE 
Q[z,E(z)} x — = 

with Q[z, E(z)\ an algebraic function of redshit z, the 
dimensionless Hubble parameter E(z) and the constant 
parameter of the model. Since we know that the uni- 
verse is expanding, dE/dz > so that we must solve 
Q[z, E(z)] = 0. Rearranging the different terms gives: 



W + y 2 ~V Ml - In (1 + rjy)] + [fi M + fi r (l + «)] 
x(l + z) 3 }y + Kln{l + riy) = [Q M + fi r (l + z)]{\ + z) 3 } 

(24) 

where we have set y — E(z) and defined : 



3 a 



(25) 



KT) 



KT] 



Kffy 



vV l + rfy 2 V 1 + v 2 y 2 

2 [Q M + fir(l + z)] (1 + zf + K ln (py + ^l + ^V) 

y 2 

with y, rj and k defined as above. The relation between 
rj and the other model parameters may be obtained as 
before and turns out to be : 



In ( i] + \/\ + rj 2 
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(29) 



which has still to be solved numerically. Note that, as 
for the simple case, also the standard MONDian vector 
model is characterized by only two quantities, namely the 
total matter density parameter VLm and the present day 
scaled Hubble constant h. 



III. TESTING THE LOW REDSHIFT REGIME 
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with h = i/o/(100/km/s/Mpc). Note that, in the rhs 
of Eq. (|2"6")l , we have used Eq. ([!§)) and set a = 1.2 x 
10~ 10 m/s 2 in agreement with the estimates coming from 
the MOND fit to the galaxy rotation curves. 

In order to reduce the number of parameters of the 
model, we can insert Eq. (|22|) into (|12| and evaluate it at 
z = 0. Remembering that E(z = 0) = 1 by definition, 
we then get : 
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Any model that aims to describe the evolution of the 
universe must be able to reproduce what is indeed ob- 
served. This is particularly true for the models we are 
considering since we are introducing a single vector field 
to get rid of both dark matter and dark energy. Match- 
ing the model with observations is also a powerful tool to 
constrain its parameters and, as an interesting byprod- 
uct, allows us to estimate some quantities common to 
every model (such as the age of the universe to and the 
transition redshift zt) to previous literature values. 

As a first test, we are here interested in exploring the 
behaviour of our models in the low redshift regime so 
that we start by only considering data probing up to 
z ~ 1.5. To this end, we therefore maximize the following 
likelihood function : 
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(30) 



where p denotes the set of model parameters. Before 
discussing in detail the term related to the SNela data, 
we concentrate on the two Gaussian priors. The first one 
takes into account the constraints on the the physical 
baryon density oj^ — flbh 2 with 



oh 8 



±Suj b = 0.0228 ± 0.0055 



as estimated in [30j and in agreement with what is in- 
ferred from the abundance of ligth elements [3l[. The 
HST Key Project [HJ has estimated the Hubble constant 
Ho using a well calibrated set of local distance scale es- 
timators thus ending up with 

h HST ±5h = 0.720 ±0.008 
as a final model independent constraint. 

A. The SNela data 

Since the first announcement [2| of the evidence of cos- 
mic speed up, the importance of the SNela Hubble dia- 
gram as a probe of the universe background evolution 
has always been clear. It has therefore become a sort of 
ground zero for every proposed cosmological model to fit 
the SNela data. To this end, one relics on the predicted 
distance modulus 2 : 



Hth{z,P) = 25 + Slog 



— (l + z)r(z,p) 

#0 



(31) 



with r(z) the dimensionless comoving distance : 



r(z,p) 



dz' 



io E(z',p) ■ 
The likelihood function is then defined as 
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where NsNeia is the total number of SNela used, A/i is a 
N SNela - dimensional vector with the values of fi b s (zi) — 
fith(zi) and C SNela is the NsNeia x N SNela covariance 
matrix of the SNela data. Note that, if we neglect the 
correlation induced by systematic errors, CsNeia is a di- 
agonal matrix so that Eq. (|3"3")l simplifies to : 



CsNeIa{p) OC exp ['XSNela (p)/ 2 ] 



(34) 



with 



XsNel 



,(p)= E 



fJ>obs(Zi) - l-Mh(zi) 
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(35) 



with Gi the error on the observed distance modulus 
l^obs(zi) for the i- th object at redshift z^. As input data, 
we use the Union SNela sample assembled in [33[ by re- 
analysing with the same pipeline both the recent SNela 
SNLS @ and ESSENCE [| samples and older nearby 
and high redshift [4j datasets. 



B. How many model parameters ? 



We use In a; and log a: to denote the logarithm base e and 10. 



As shown by Eqs. ([24|) and ([28]) . in order to deter- 
mine the cosmic dynamics, one has to set the value of 
the present day total matter density parameter f2jvf- On 
the other hand, the Gaussian prior on u)b in the likeli- 
hood function (|30p and of the distance priors introduced 
later asks for discriminating between the baryons only 
and total matter physical densities. It is therefore worth 
wondering how Qm and Qb are related. Since our the- 
ory reduces to MOND in the low energy limit and MOND 
does not need any cold dark matter on the galactic scales 
(hence no other matter than the visible one), we could ar- 
gue that Qjvf — ilb should hold. On the other hand, large 
amounts of missing matter are needed in order MOND to 
reproduce the observed phenomenology on cluster scales 
[34 ] . It was therefore postulated [35[ that massive neu- 
trinos (with m v ~ 2 eV) can play the role of dark mass 
in galaxy clusters. Moreover, solar and atmospheric neu- 
trinos experiments [36| have shown that the three active 
neutrinos from the standard model of particle physics 
mix their flavours which is only possible if they are mas- 
sive. Finally, it is also worth noting that it has been 
claimed (37j that a single sterile neutrino with mass in 
the range 4-6 eV is better suited to explain the results 
of the Miniboone experiment, while a 11 eV sterile neu- 
trino has been indeed advocated [38| in order to solve the 
MOND problems on cluster and cosmological scales. 

Massive neutrinos decoupled at ~ 1 MeV and since 
last scattering they have been non relativistic particles 
so that, from a cosmological point of view, they behave 
exactly as matter. Assuming three families of degenerate 
neutrinos, the total matter density parameter will there- 
fore read : 
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Par x (x) Xmed 68% CL 



95% CL 



fib 0.13 0.12 (0.05, 0.23) (0.02, 0.29) 

logm„ 0.24 0.40 (-0.16,0.55) (-0.88,0.63) 

h 0.700 0.700 (0.694, 0.706) (0.689, 0.711) 

fit 0.15 0.13 (0.07, 0.26) (0.02, 0.31) 

log m„ -0.09 0.35 (-1.10,0.53) (-2.40,0.64) 

h 0.700 0.700 (0.694, 0.706) (0.688, 0.712) 



TABLE I: Summary of the results of the likelihood analysis 
including SNela and Gaussian priors on Wb and h. The upper 
and lower part of the table refers to the simple and standard 
MONDian vector model respectively. 



M 



fi h 



3m„ 
94h 2 eV 



(36) 



so that the number of model parameters is increased by 
one updating from (fijif,h) to (Q;,, m„, h). 

It is worth stressing, however, that discriminating be- 
tween Qm as a single quantity and Qm as function of 
(fib, m u , h) is only possible if the data at hand depend 
on them separetely. To understand this point, let us 
consider the SNela Hubble diagram. In order to fit this 
dataset, we just need the dimensionless Hubble parame- 
ter E(z) which is obtained by solving, e.g., Eq. (|2"4")) . To 
this end, we just have to set the value of f2 M so that all 
the SNela likelihood will be a function of (S1m,/i) only 
and we cannot set any constraint on (flb,m v ). That is 
why we have added the Gaussian priors on ujb and h 
in order to make the likelihood explicitly dependent on 
fib- However, since the likelihood is mainly driven by 
the SNela term, one can forecast that the constraints on 
(fib, rn v ) will be quite weak because of the degeneracy 
being only partially broken. 



C. Results 

In order to maximize the likelihood function (I30|) . we 
run a Markov Chain Monte Carlo (MCMC) algorithm to 
efficiently explore the parameters 3 space. We use a sin- 
gle chain with ~ 100000 points reduced to ~ 3100 after 
burn in cut and thinning. The histograms of the values of 
each parameter are then used to infer the median and the 
mean and the 68 and 95% confidence ranges summarized 
in Table I. The best fit parameters (i.e., the set maximiz- 
ing the likelihood) for the simple MONDian vector model 
turns out to be : 

(fl h ,\ogm v ,h) = (0.056,0.54,0.70) , 



while the standard case gives : 

(fl b ,logm v ,h) = (0.050,0.55,0.70) . 

For both models we get XsNeia/d-°-f- — 1-02 thus in- 
dicating a very good agreement. Moreover, the physical 
baryon density reads u>b = 0.0276 (0.0246) for the simple 
(standard) model in good agreement with the observed 
value well within the ler error, while both models predict 
values for h very close to the HST Key Project result. 
It is worth noting that a baryon only universe may be 
safely excluded since the reduced x 2 values are of order 
2.15 for both cases thus definitevely ruling out models 
with no massive netrinos. 

A look at Table I shows that the constraints on both 
fib and logWjy are quite weak. Moreover, the best fit 
values of both parameters radically differ from their me- 
dian values. This is, however, not an unexpected result. 
As we have said, should the prior on u>b be neglected, the 
model should collapse into a two parameter one with VLm 
replacing the set (fib, logm„). The prior on Ub thus helps 
in discriminating between the two, but, from the point 
of view of maximizing the likelihood, it is of little help. 
Indeed, the only reliable constraint is on flm reading : 



(fl M ) = 0.28 , fl 



M,med 



= 0.28 



3 Note that, hereafter, we use logm„ rather than m v as neutrino 
mass parameter since this choice allows to explore a wider range. 



68% CL : (0.25,0.31) , 

95% CL : (0.23,0.34) , 
for the simple MONDian vector model and 
(Am) = 0.28 , flM,med = 0.28 , 

68% CL : (0.25,0.31) , 

95% CL : (0.22,0.33) , 

for the standard case. Note that these values are in 
very good agreement with typical estimates from previ- 
ous analyses of comparable datasets [2, LLCj, [33J . It is also 
worth noting that the results are almost fully indepen- 
dent on the functional expression adopted for T(1C) which 
is an expected consequence of the two models matching 
each other in order to fit the same SNela Hubble diagram. 
Provided that the above constraints on Qm are met, we 
can choose whatever value of fl b (for a given h) and then 
find a corresponding log m v value giving rise to a model 
with a given likelihood value, i.e. C depends only on 
(Om, h). As a consequence, it is therefore not surprising 
that the constraints on (fib, logm„) are so weak. 

With this caveat in mind, it is nevertheless interesting 
to look at the neutrino mass. Converting the constraints 
on \ogm v into constraints on m„ (in eV), we get : 
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(m„) = 2.3 , m v , m ed = 2.5 , 

68% CL : (0.7,3.6) , 

95% CL : (0.1,4.2) , 
for the simple MONDian vector model and 
(m„) = 2.0 , m u>me d = 2.2 , 

68% CL : (0.1,3.4) , 

95% CL : (0.0,4.4) , 

for the standard case. As yet stated, atmospheric and so- 
lar neutrino experiments have shown that the three fami- 
lies of standard model neutrinos are massive, but they are 
unable to put any constraints on their exact masses be- 
ing only sensitive to mass squared differences. An upper 
limit on the mass may instead be set from the study of the 
tritium j3 decay. By this method, the Mainz - Troitz ex- 
periment [39| was able to find m v < 2.2 eV. The median 
m v values quoted above are smaller than this upper limit, 
while the 68 and 95% confidence ranges do indeed suggest 
that it is possible to fit the data equally well with still 
lighter neutrinos. On the other hand, it is worth stressing 
that such estimates rely on our assumption that three de- 
generate massive neutrinos are present so that their total 
density parameter reads fl v = 3m„/94ft, 2 . It has, how- 
ever, been claimed [33] that a single sterile neutrino with 
mass in the range 4-6 eV is better suited to explain the 
results of the Miniboone experiment. Should this be the 
case, the constraints above should be multiplied by three 
thus giving a single sterile neutrino with mass in the 95% 
confidence range (0.3, 12.6) eV. It is worth noting that a 
11 eV sterile neutrino has been indeed advocated [38| in 
order to solve the MOND problems on cluster and cos- 
mological scales. It is therefore tempting to investigate 
whether the results in [38[ still hold for our vector models 
since on cluster scales both theories reduce to MOND. 

The MCMC algorithm makes also possible to infer con- 
straints on some interesting derived quantities. To this 
end, one just has to evaluate a given function /(p) for 
each point of the chain and then use the values thus ob- 
tained as yet done for the parameters p. Table II sum- 
marizes the results for the present day deceleration pa- 
rameter go j the transition redshift zt obatined by solving 
q(zx) = and the age of the universe to estimated as 

to = tH [ 77—. — ttttt (37) 
Jo (l + z)E(z) 

with t H = 9.78/i _1 Gyr the Hubble time. 



Par x (x) x med 68% CL 95% CL 

go -0.57 -0.58 (-0.62, -0.53) (-0.65, -0.49) 

z T 0.73 0.72 (0.65, 0.80) (0.58, 0.88) 

t 13.71 13.70 (13.43, 14.00) (13.18, 14.31) 

go -0.58 -0.58 (-0.62, -0.53) (-0.65, -0.49) 

zt 0.73 0.72 (0.64, 0.81) (0.57, 0.88) 

t 13.72 13.71 (13.42, 14.02) (13.16, 14.29) 

TABLE II: Constraints on derived quantities (to in Gyr) from 
the chains obtained fitting the SNela with Gaussian priors on 
u>t and h. The upper and lower part of the table refers to the 
simple and standard MONDian vector model respectively. 



As a first issue, let us consider the value of go- Its esti- 
mate is typically model dependent since it comes out as 
a derived quantity given a model parametrization. In 
order to escape this problem, one may resort to cos- 
mographic analyses based only on Taylor expanding the 
scale factor. Using this approach, Cattoen and Visser 
[iol ] have found values between g = —0.48 ±0.17 and 
go = —0.75 ±0.17 depending on the details of the method 
used to fit the SNLS dataset. A similar analysis but us- 
ing the GRBs as distance indicators allowed Capozziello 
and Izzo [4l[ to find values between go = —0.94 ± 0.30 
and go = —0.39 ± 0.11 still in accordance with our esti- 
mates. A different approach has been instead adopted by 
Elgar0y and Multamaki [42j advocating a model indepen- 
dent parametrization of q(z). Depending on the SNela 
sample used and the parametrization adopted, their best 
fit values for go range between -0.29 and -1.1 in good 
agreement with the estimates in Table II. 

On the contrary, there is some conflict with the previ- 
ous estimates of the transition redshift Zt- For instance, 
using the Gold SNela sample and linearly expanding q(z), 
Riess et al. [J] found zt — 0.46 ± 0.13 in agreement with 
the updated result zt — 0.49lo'o7 obtained by Cuhna 
[43l ] using the Union sample. Although there is a possi- 
ble marginal agreement within the 95% confidence range, 
we nevertheless consider this unsatisfactory result not a 
serious flaw of our models since the estimate of zt is 
strongly model dependent so that it is not possible to 
decide whether the disageement is with the data or with 
the fiducial model used to fit the data. 

Finally, we note that the age of the universe is in 
agreement with previous estimates in the literature. Fit- 
ting the WMAP5 data and the SNLS SNela sample 
with a prior on the acoustic peak parameter gives [10( 
to = 13.73 ± 0.12 Gyr in almost perfect agreement with 
the results in Table II. Moreover, the shape of the colour - 
magnitude diagram of globular clusters provides a model 
independent estimate, namely to — 12.61^'g Gyr [44| . 
still in considerable good agreement with our values. 

Summarizing, the very good fits to the SNela data 
and the agreement between observed and predicted de- 
rived quantities make us confident that both the simple 
and standard MONDian vector models successfully re- 
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produces the data thus being viable alternatives to the 
usual dark energy models in the low redshift regime. 



IV. PROBING THE HIGH REDSHIFT REGIME 

The SNela Hubble diagram and the Gaussian priors 
on uib and h allow us to test the behaviour of the MON- 
Dian vector models only over the redshift range probed 
by these data. Considering that the farthest SN has red- 
shift z ~ 1.6, it is worth wondering whether the mod- 
els work well for higher z. To answer this question, we 
change the likelihood function adding a further term : 



£(p) oc Aow(p) x £d P (p) 



(38) 



with £iow(p) the likelihood term related to low redshift 
data given by Eq. (|30|) . while £d P (p) is the distance priors 
dependent term which we detail below. 



A. The distance priors 

While SNela probe only the background evolution of 
the universe over a limited redshift range (up to z ~ 1.5), 
the CMBR anisotropy spectrum and the matter power 
spectrum measured by galaxy surveys make it possible 
to both constrain the dynamics of perturbations and test 
the model up to the last scattering surface. However, 
such a program is both theoretically difficult and com- 
putationally demanding. On the one hand, we must de- 
velop the full theory of perturbations to compute the 
Ci coefficients of the multipole expansion of the CMBR 
anisotropies and the matter power spectrum P(k). While 
this is yet done for the standard dark energy models, 
a full perturbation theory for the vector models we are 
considering is still to be developed. On the other hand, 
even for popular dark energy models, computing both the 
Ci and P(k) represents a bottleneck for the algorithms 
matching data with theories. Fortunately, many features 
of the CMB and matter power spectra may be expressed 
as a function of a limited subset of quantities which are 
instead easy and straightforward to compute. Motivated 
by this consideration, it has become popular to summa- 
rize the main constraints coming from CMBR and matter 
power spectra in what are defined distance priors . 

In the analysis of the WMAP5 data [lo| . Komatsu et 
al. demonstrated that most of the information in the 
WMAP power spectrum may be summarized in a set of 
constraints on the following quantities : 



(i.) the physical baryon density: 



uj b = fl b h 2 (39) 

with fib the density parameter of baryons only; 

(ii.) the redshift z^s to the last scattering surface that 
we approximate as [45| : 

z LS = 1048(l+0.00124 W -°- 738 )[l+ 5l (^) W § ( ^ ) ] , (40) 



. . 0.0738wr 0238 



92M 



1 + 39.5cf 763 ' 



0.560 
l + 21.1wi- 81 



(41) 



(42) 



having denoted with lo^i = Clj^h the total matter 
physical density; 

(iii.) the acoustic scale [46| : 



Ia 



ir(c/H a )r(z LS ) 
r s (z L s) 



(43) 



with r s (zLs) the size of the sound horizon at the 
decoupling epoch given by : 



r s {a LS ) 



c/H 



da 



Jq a 2 E(a)\ / l + Ra 
being ols = (1 + zls)^ 1 and R = 3^/4^; 
(iv.) the shift parameter defined as [46| : 

TZ = \JVL M r{z LS ) . 



(44) 



(45) 



While this set of constraints relies on the CMBR data, 
none of them has to do with the matter power spectrum. 
To overcome this problem, Eisenstein et al. [471 ] intro- 
duced the acoustic peak parameter defined as 



A= v^m Dv{Zm) 



CZ 7 -, 



(46) 



where z m is the median redshift of the galaxy survey used 
to extract the matter power spectrum and the volume 
distance is given by : 



4 Actually, not all theese quantities are indeed distances. Never- 
theless, the set of constraints is collectively referred to with this 
name since most of them are easily related to a distance. 



Dv(z) = 



H(z) 



cr 2 (z) 



1/3 



(47) 
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Par x (x) x med 68% CL 95% CL 

fi b 0.0394 0.0394 (0.0385, 0.0404) (0.0374, 0.0413) 

log m„ -0.036 -0.037 (-0.048,-0.025) (-0.059,-0.012) 

h 0.753 0.753 (0.749, 0.757) (0.745, 0.761) 

fi b 0.0384 0.0384 (0.0372, 0.0394) (0.0363, 0.0402) 

logm„ 0.005 0.006 (-0.016,0.025) (-0.031,0.041) 

h 0.753 0.753 (0.749, 0.758) (0.745, 0.761) 



TABLE III: Summary of the results of the likelihood analysis. 
The upper and lower part of the table refers to the simple and 
standard MONDian vector model respectively. 

The analysis of the correlation function of a sample of 
more than 46000 LRGs from the SDSS allowed Eisein- 
stein et al. to estimate : 

A= 0.469 ±0.017 

thus offering another valuable constraint on the model 
parameters 5 . However, a more recent analysis of the sev- 
enth SDSS data release allowed Sanchez et al. [30j to 
measure the correlation function along both the radial 
and tangential directions thus allowing a more detailed 
treatment. Adding the CMBR data in a joint analysis, 
the authors then provided a new set of distance priors 
adding to the four quantities quoted above a new one 
defined as 

G(z m ) = r(z M ) x [E(z m )] Q - s . (48) 

In order to take into account the distance priors, we in- 
troduce the following likelihood function : 

r t \ 1 ( Ad P C dp A d P \ r . n , 
**CP) - (27r)6 / 2 | C ^|i/2 >< «P [ 2 ) (49) 

where A,j p is a five dimensional vector whose i - th ele- 
ment is given by A^j = Xi t0bs — %i,th with the subscript 
obs and th denoting the observed and theoretically pre- 
dicted values. The label i runs from 1 to 5 referring, 
respectively, to u b , zls, I a, ft and G(z m ). Finally, we 
follow Appendix B of [30] to set the observed values and 
the covariance matrix of the distance priors. 



B. Fitting SNela and distance priors 

We now run our MCMC algorithm to explore the mod- 
els parameters space now maximimizing the likelihood 
(f58"|) taking care of both the SNela and distance priors 



5 Note that A does not depend on h. 



datasets. The median and mean values and 68 and 95% 
confindence ranges are summarized in Table III for both 
the simple and standard MONDian vector model. The 
best fit parameters, i.e. the values of (Oj,, logm„, h) that 
maximizes the likelihood, are not forced to be equal to 
the median values because of correlations among the pa- 
rameters. Indeed, we find : 

(fi 6 , logm„, h) = (0.0392, -0.038, 0.753) 
for the simple MONDian vector model, while it is : 

(Q b , logm„, h) = (0.0387, -0.001, 0.753) 

for the standard case. Both models, however, do not 
offer a good performace in fitting the data. Indeed, for 
the best fit parameters, we get : 

X% N eiJd-o.f = 1.32 , oj b = 0.0222 , z LS = 1081.7 , 

Z A = 300.0 , ft =1.52 , G(z m ) = 1424 , 
for the simple case and 

XsNeiJd-o-f = 1-31 , uj b = 0.0219 , z LS = 1082.3 , 

l A = 300.4 , ft =1.54 , G(> m ) = l424 , 

for the standard model. It is immediately clear from the 
unusually high reduced XsNeia va lue signals that some- 
thing is going wrong with fitting the SNela Hubble dia- 
gram. Indeed, with d.o.f. = NsNeia—3 — 304, a reduced 
XsNeia/d-O-f- — l-3hasjust a tiny ~ 3xl0 -5 probability 
to occur so that we can safely conclude that both mod- 
els are not correctly fitting the SNela Hubble diagram 
data. This conclusion is further enforced comparing the 
above best fit distance priors whose median values and 
standard deviation 6 are as follows : 

Lu b = 0.0228±0.0055 , z LS = 1090.1±0.9 , l A = 301.6±0.7 , 

ft = 1.701 ±0.018 , G{z m ) = 1175 ±21 . 

While the values of (uj b , z^s, I a) are in reasonable agree- 
ment, there are strong discrepancies for both the shift ft 
and the G(z m ) parameters. Considering the 68 and 95% 
confidence ranges in Table IV does not ameliorate the 
comparison so that we must conclude that the model is 
unable to fit both the SNela and distance priors dataset. 



6 Actually, the standard deviation is not a good estimator of the 
uncertainty on the distance priors because it does not take into 
account the correlations among them. However, it gives an idea 
of the discrepancy between predicted and observed values. 
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Par x (x) x m ed 68% CL 95% CL 

ui b 0.0224 0.0223 (0.0219, 0.0229) (0.0213, 0.0234) 

z LS 1081.6 1081.6 (1081.1, 1082.0) (1080.7, 1082.6) 
l A 299.9 299.9 (299.2, 300.5) (298.6, 301.2) 
Tl 1.524 1.523 (1.519, 1.529) (1.514, 1.533) 

G{z m ) 1424 1424 (1416, 1431) (1407, 1439) 

u b 0.0218 0.0218 (0.0211, 0.0223) (0.0206, 0.0228) 

z LS 1082.6 1082.5 (1081.8, 1083.4) (1081.2, 1084.0) 
Ia 300.6 300.6 (299.9, 301.4) (299.1, 302.2) 
TZ 1.545 1.545 (1.543, 1.549) (1.539, 1.551) 

G(z m ) 1424 1424 (1416, 1432) (1409, 1439) 

TABLE IV: Constraints on the predicted distance priors pa- 
rameters. The upper and lower part of the table refers to the 
simple and standard MONDian vector model respectively. 

It is worth investigating why this happens. Actually, 
a hint is given by noticing that the most discrepant pa- 
rameters are those involving E(z) both directly, as for 
G{z m ), or indirectly through an integral, as both 1Z and 
XsNeia- Taking the ACDM model as a comparison, we 
indeed find that both the simple and standard MONDian 
vector models systematically underestimates E(z). As a 
consequence, the dimensionless comoving distance r(z) 
turns out to be overestimated thus leading to /i(z) be- 
coming increasingly higher than the concordance model 
prediction as the redshift increases in accordance with the 
result that /i(z) is larger than /z {, s for SNela with z > 1. 
This also explains why G(z m ) gets larger than observed, 
while the situation is different with the shift parameter. 
Indeed, according to Ea.(|35]). a larger r(z) should lead 
to a larger TZ, while we observe the opposite result. This 
can, however, be easily explained considering that TZ is 
also proportional to J7m and we get : 

(Sl M ) = 0.091 , n Mtmed = 0.091 , 

68% CL : (0.089,0.094) , 

95% CL : (0.086,0.096) , 
for the simple MONDian vector model and 

(Sl M ) = 0.095 , n M ,med = 0.095 , 

68% CL : (0.094,0.097) , 

95% CL : (0.092,0.098) , 

for the standard case. These values are much smaller 
than the typical Am — 0.26 @,[HJ[III values obtained in 
literature thus explaining why our shift parameter turns 
out to be so small notwithstanding the higher r(zLs)- 



Par x (x) x med 68% CL 95% CL 





-0.861 


-0.862 (-0.865, -0.858) (-0.869, -0.855) 


Zt 


1.71 


1.71 (1.68, 1.73) (1.66, 1.76) 


to 


16.93 


16.93 (16.84, 17.01) (16.76, 17.10) 


qa 


-0.856 


-0.856 (-0.859, -0.854) (-0.861, -0.852) 


ZT 


1.67 


1.67 (1.65, 1.68) (1.64, 1.70) 


to 


16.77 


16.77 (16.66, 16.87) (16.60, 16.98) 



TABLE V: Constraints on derived quantities (with to in Gyr) . 
The upper and lower part of the table refers to the simple and 
standard MONDian vector model respectively. 

For completeness, we also summarize in Table V the 
values of the present day deceleration parameter qo, the 
transition redshift zt and the age of the universe to. 
These results strengthen our conclusion that both mod- 
els, while performining well in the low redshift regime, are 
actually quite poor in reproducing data probing higher z. 
Indeed, the values of go are in reasonable agreement with 
the results quoted in Sect. IIIC, even if our qo values are 
somewhat extreme. On the contrary, there is a clear dis- 
agreement with the previous estimates of the transition 
redshift zt with our results being very high (also outside 
the range probed by SNela Hubble diagram). Such a re- 
sult furtherly signals that indeed the expansion rate and 
hence the transition from acceleration to deceleration of 
our models is too slow. Another evidence in favour of 
this interpretation is provided by to which turns out to 
be in strong disagreement with both the WMAP5 and 
globular clusters estimate. Indeeed, Eq. ([37|) shows that, 
should we underestimate E(z), the age of the universe 
turns out to be higher which is just what happens for 
both MONDian vector models. 



C. A problem with the models or the data ? 

The large reduced XsNeia values, the disagreement be- 
tween the best fit predicted and observed distance priors 
and the unacceptably high zt and to results should make 
us conclude that both MONDian vector models are un- 
able to fit both the low and high redshift data. Com- 
paring the results in Table I and III, however, makes it 
evident that the two fitting procedures select very differ- 
ent regions of the parameter space. It is therefore worth 
wondering whether there is a problem with the distance 
priors data rather than with the models. 

To this end, it is worth stressing that, alth oug h widely 
used in the recent literature (see, e.g., ffl. Hoi. l33l liH. |49]| ) . 
their estimate is actually model dependent. Indeed, in or- 
der to obtain their central values and covariance matrix, 
one first fits a given model (typically, the concordance 
ACDM one) to the full CMBR anisotropy and galaxy 
power spectra dataset using a Markov Chain Monte Carlo 
(MCMC) method to sample the posterior probability. 
This same model is then used to compute the distance 
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FIG. 1: Effective EoS (left) and dark energy density parame- 
ter (right) for the standard MONDian vector model. We set 
(Qt,h) — (0.04,0.70) and consider three values for logm„, 
namely -1.0 (short dashed), 0.0 (solid) and 0.5 (long dashed). 



priors along the chain and then the sample thus obtained 
is analysed to infer the covariance matrix. As correctly 
stressed in 50], this procedure relies on three main as- 
sumptions : 

1. the posterior probability from the CMBR 
anisotropics and galaxy power spectra is cor- 
rectly described by the distance priors, i.e. no 
information is lost when giving away the full 
dataset in favour of the simplified one; 

2. the mean and covariance matrix of the distance pri- 
ors parameters do not change when the model space 
is enlarged, i.e. the choice of the fiducial model does 
not impact the estimate of the priors; 

3. these summary parameters are weakly correlated 
with the input fiducial model parameters so that 
they can be used as independent constraints. 

While the first and third assumptions have indeed been 
verified [H, H3] , the second hypothesis can not be fully 
addressed. On the one hand, one could argue that, un- 
less a significative amount of dark energy in the early 
universe is present, all reliable models should match at 
high redshift, when the effect of dark energy fades away. 
As a consequence, using the ACDM as a fiducial model 
to compute quantities that mainly depend on the high 
redshift behaviour should not affect the final estimate. 
However, even if the background evolution could be the 
same, it is still possible that the dynamics of perturbation 
is radically different as is the case with modified gravity 
theories. This is indeed what also happens for the MON- 
Dian vector models we are considering. Fixing the model 
parameters to the best fit values from the fit to the low 
redshift data only, we can easily compare the luminosity 
distance with that for a fiducial ACDM scenario. As ex- 
pected (since they both reproduce the same data), the 
two can hardly be discriminated over the range probed, 
while the Hubble parameter only differ by a few percent. 
While the background evolution is therefore comparable, 
perturbations evolve in a completely different way be- 
cause of the presence of the vector field. As such, one can 
not be sure that the information contained in the CMBR 
spectrum may still be summarized in the distance priors 
quantities as for the ACDM model. A similar discussion 
also applies for the the prior on G(z m ). Moreover, it is 



FIG. 2: Same as Fig.[TJ but for the simple MONDian vector 
model. The small ripples are due only to numerical errors. 



worth stressing that this latter quantity actually depends 
on data probing an intermediate redshift range so that 
one can not rely anymore on the fading of dark energy 
at high z as is typically done for quantities as, e.g., the 
acoustic scale I a and the shift parameter 1Z. 

Investigating the impact of these problems on the es- 
timate of the distance priors is outside our scope here. 
We nevertheless stress that, because of the above consid- 
erations, one can not safely reject the MONDian vector 
models because of the disagreement with the distance pri- 
ors values. As a conservative conclusion, we are therefore 
forced to only report the results being unable to decide 
whether the problem is with the data or with the fiducial 
model used to retrieve them. 



D. The effective EoS and the high z limit 

An alternative way to compare the proposed MON- 
Dian vector models with standard dark energy models in 
both the low and high redshift regime may be obtained by 
considering the effective EoS. Indeed, from the point of 
view of the background evolution, our models are equiv- 
alent to a cosmological scenario made out of dust matter 
and a dark energy with an effective EoS given by : 



W eff(z) 



2 d\riE(z) n M (l + z) 
3dln(l + z) E 2 (z) 
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1 - 



n M {i + zf 



E 2 (z) 



(50) 



so that the dark energy density parameter reads : 



^de{z) 



i - n 



E 2 (z) 



- exp 



eff(z') 



1 



dz 



(51) 



Figs.Q]and[2]show w eff(z) and SIde{z) for the standard 
and simple MONDian vector models setting the baryon 
density parameter to the fiducial value Qb — 0.04 and 
varying the neutrino mass m„. Note that the value of 
Am in Eqs.([5"0|) and (fSTj) is related to Sib and m v through 
Eci. (f3"S|) so that changing m v is the same as varying fljj. 

Both figures shows that the value of 77i„ set the present 
day EoS with w e f t(z = 0) increasing as a function of m„, 
even if there is a saturation at small m v as can be easily 
understood considering that the smaller is m v , the less 
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is the impact of massive neutrinos in the energy budget. 
It is worth noting that for logm„ ~ 0.5, it is w e ff(z — 
0) ~ —1.0 and dw e ff/dz takes very small values over 
the redshift range (0, 1), i.e. we recover (for both cases) 
a present day cosmological constant. It is therefore not 
surprising that the fit to the SNela dataset points towards 
such values of log m„ since these make both MONDian 
vector models as similar as possible to the ACDM one 
over the range probed by SNela themselves. 

On the contrary, in the high redshift regime, the neu- 
trino mass only plays a marginal role in determining the 
value of w e ff(z) which stays almost constant at value 
close to the dust one w e ff = 0, i.e. the effective EoS 
approximately tracks the matter term. However, the 
amount of dark energy is different depending on m„ with 
lower values of the neutrino density giving rise to a larger 
SIj)b(z) in the high z regime. Such a behaviour is some- 
what counterintuitive since one expects that the effective 
dark energy fades away in the early universe in order to 
recover the matter dominated epoch. Actually, one must 
take into account that our matter term is made out of 
baryons and massive neutrinos only. The lower is the 
neutrino mass, the higher must be the contribute of the 
effective dark energy (that in this regime behaves as mat- 
ter being w e ff ~ 0) to compensate for the missing cold 
dark matter. Such a result can also be forecasted going 
back to Eq.© and noting that, for our best fit models, 
we find T) » 1 so that, for z >> 1, it is K >> 1 too. 
Let us then consider the simple MONDian vector model. 
Inserting Eq. ([22| and taking the limit for r]E(z) >> 1, 
we approximately get : 



n M (l + zf ~E 2 + bxrjE ~ E 2 



so that we indeed recover the usual expression for the 
Hubble parameter for a matter dominated universe with 
Ojf the total (baryons and massive neutrinos) matter 
term (having neglected radiation). From the point of 
view of the effective dark energy formalism, if we reduce 
the neutrinos contribution, we must add a further term 
acting as matter and this is indeed provide by the effec- 
tive dark energy. This also explains why Qde(z >> 1) 
increases decreasing log m v . A similar discussion also ap- 
plies to the standard MONDian vector models so that in 
both cases the final scenario is the one of a universe where 
the matter term disappears to be replaced by a dark en- 
ergy fluid acting approximately as matter. This explains 
why the inclusion of the distance priors (probing the high 
z regime) pushes the best fit towards smaller log m„ since, 
in this case, both MONDian models recovers the usual 
Fricdmann models which are known to successfully fit 
these high redshift probes. As yet said, however, de- 
creasing logm„ makes w e ff(z = 0) ^ — 1 thus worsening 
the fit to the SNela Hubble diagram. 



V. CONCLUSIONS 

The astonishing successes of MOND on the galactic 
scales and the emergence of a relativistic theory playing 
the role of its counterpart on cosmological scales have re- 
newed the interest in the search for a possible common 
explanation of both dark matter and dark energy phe- 
nomenology. Vector theories are another way of recov- 
ering MOND in the low energy limit so that it is worth 
wondering whether they can also offer an elegant way 
of speeding up the cosmic evolution without the need of 
any dark energy source. The ignorance about the form 
of the vector field Lagrangian may be bypassed relying 
on the link between the function J-(K.) and the MOND 
interpolating function /x(a/ao). Since dark matter is no 
more present, one should postulate the presence of mas- 
sive neutrinos in order to fill the gap between the total 
matter density parameter and the baryons density alone. 
Moreover, such massive neutrinos are also advocated in 
order to reconcile the results of solar and atmospheric 
neutrino experiments on the flavour mixing with the pre- 
dictions of the standard model of particle physics. 

Motivated by these considerations, we have there- 
fore investigated the viability of two different MONDian 
vector models characterized by !F(JC) expressions corre- 
sponding to the simple and standard MOND interpo- 
lating function. To this aim, we have first fitted them 
against the Union SNela Hubble diagram using Gaussian 
priors on the physical baryon density u>b and the present 
day Hubble constant h in order to break the (fib, \ogm u ) 
degeneracy. Both models performs quite well giving a 
perfect agreement with the SNela data and previous es- 
timates of the total matter density parameter fijvf, the 
deceleration parameter qq and the age of the universe to- 
Moreover, the (weak) constraints on the neutrino mass 
are consistent with the upper limits set by the Mainz - 
Troitz experiment assuming three families of degenerate 
neutrinos. Should, instead, a single sterile neutrino be 
the mass dominant component, its estimated mass is in 
agreement with the constraints from the Miniboone ex- 
periment also falling in the right range advocated to solve 
MOND problems on cluster scales. 

In order to investigate the high redshift behaviour of 
the models, we have repeated the likelihood analysis 
adding the extended set of distance priors. It turns out 
that both models should be rejected since they provide 
now a poor fit to the SNela data and strongly disagree 
with the observed shift 1Z and G(z m ) parameters. Dis- 
tance priors are, however, estimated through a model 
dependent procedure so that one can not safely rely on 
them to exclude models that are radically different from 
the fiducial one used to extract the constraints on 1Z and 
G(z m ). As a consequence, we are unable to conclude 
whether the disagreement is a failure of the MONDian 
vector models or an expected outcome of the different 
evolution of perturbations with respect to usual dark en- 
ergy models because of the action of the vector field. 

It is worth remembering that the search for a cosmolog- 
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FIG. 3: Left : the g(x) function for the best fit standard MON- 
Dian vector model (solid line) compared to the Lue & Stark- 
man proposals. Right : relative deviation of the proposal g(x) 
from that of our model. Short and long dashed lines refer to 
g(x) given by Eqs. (f52)l and (f53|) respectively. 



ical counterpart to MOND has a long history with vector 
models being only the most recent proposal. Applying a 
procedure which is simply the generalization of the clas- 
sical derivation of the Friedmann equation from the New- 
tonian force law, but now starting from the MOND force 
law, Lue & Starkman [5ll | derived an expression for the 
Hubble parameter without referring to any underlying 
modified gravity theory. Defining g(x) = H 2 /Hq and 
x = + z ) 3 , they indeed find : 



(Plnx + c 2 )a; 2 / 3 x < x 



(52) 



X + C\X 



2/3 



with ci = C2 + 3/3[ln (3/3) — 1]. Moroever, they also pro- 
posed a modified (by hand) version of this expression to 
better account for the low redshift behaviour : 



7(s) 



/3x 2 / 3 In (1 + z) 



x < 0.1 



0.1 < x < (3/3) 3 



x + 3/3[ln (3/3) - l]x 2 / 3 x > (3/3) 3 



(53) 

Following [511, we set (/3, x c , Q A ) = (15, 7 x 10 4 , 0.7) with 
c 2 = and show in Figs.[3] and dhow Eqs.lO and (j5"3")l 
compares to the standard and simple MONDian vector 
models. As it is better appreciated from the right pan- 
els, both Eqs. (p)2")) and (j53"|) agree reasonably well with 
the g(x) expression of our MONDian models for very 
high x, that is in the early universe. This is an ex- 
pected result since, in this redshift range, both our mod- 
els and the phenomenological Lue & Starkman proposals 
recover the typical General Relativity matter dominated 
scenario. On the contrary, for small x, i.e. in the very 
low z regime, only the modified expression (|53| matches 
reasonably well those for the standard and simple MON- 
Dian vector models as a consequence of both mimicking 
an effective cosmological constant. In the intermediate 
region, however, matching the Lue & Starkman mod- 
els to our own is impossible signalling that the strategy 
adopted by these authors is unable to trace the tran- 
sition region to General Relativity. However, it is also 



FIG. 4: Same as Fig. [3] but for the simple model. 



worth stressing that such a disagreement could be ex- 
pected since the Lue & Starkman procedure relies on the 
assumption that, whatever the underlying modified grav- 
ity theory leading to (|52|) or (|53|) is, the Birkhoff theorem 
still holds. This is not the case for vector theories ^52] so 
that the two approaches differ from the very beginning. 

As yet quoted in the introduction, much interest has 
been devoted to the TeVeS theory as a relativistic MOND 
counterpart. In particular, Skordis et al. (53| has also 
investigated the growth of structure in TeVeS also com- 
puting the CMBR anisotropy spectrum. It turned out 
that the inclusion of a cosmological constant term and of 
massive neutrinos with m v ~ 2 eV may lead to a reason- 
ably good agreement with the data, although a detailed 
fitting has not been performed. Since both the vector 
and scalar fields present in TeVeS do not contribute sig- 
nificantly to the dynamics, it is likely that this TeVeS + A 
model matches well both the SNela Hubble diagram and 
the distance priors. However, comparing the Skordis et 
al. model to our MONDian vector theories is not possible 
given the radical differences betwee the two approaches. 
Indeed, TeVeS needs the scalar field to act as a dark 
matter like term on galactic scales, while the vector field 
boosts the growth of perturbations during the A driven 
background expansion. On the contrary, in our approach, 
the vector models modify the low energy limit Poisson 
equation thus recovering the MOND - like behaviour, but 
also originates the cosmic speed up. In a sense, our ap- 
proach is more economical claiming for the lowest possible 
number of ingredients. 

The positive results obtained in the low redshift regime 
may be considered only as a ground zero level analysis. 
Indeed, fitting the SNela Hubble diagram only tells us 
that the two considered MONDian vector models predict 
the correct background evolution over the redshift range 
probed by the data, i.e. up to z ~ 1.5. Needless to 
say, more tests are needed in order to assess the viabil- 
ity of these models. On one hand, we can still investigate 
the background evolution by extending the redshift range 
through the use of the GRBs Hubble diagram [54|. On 
the other hand, a more significative and demanding test 
should be fitting both the galaxy power spectrum and 
the CMBR dataset. However, both these tasks are quite 
daunting from the theoretical point of view. Indeed, as 
far as we know, the CMBR spectrum has never been 
computed for the class of vector theories we are consid- 
ering so that one should first write down the full set of 
perturbation equations to then modify a numerical code 
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like CAMB . Somewhat easier (even if still difficult) 
is to deal with the problem of growth of structures in 
vector models. Indeed, such a study has although been 
performed in [5t| where the authors considered a simple 
power -law choice for !F(JC). Although a detailed fitting 
to the data was not performed, these authors have convic- 
ingly shown that the new degrees of freedom sourced by 
the vector field may indeed boost the growth of structure 
even in absence of any dark matter. This result is a good 
starting point for testing our proposed MONDian vec- 
tor models provided one accordingly changes the !F(K.) 
function and introduces massive neutrinos into the game. 
Note that, since we want to recover MOND on galactic 
scales, we are postulating no cold dark matter so that the 
observed galaxy power spectrum should be matched to 
the predicted one with a bias parameter determined by 
the clustering properties of the massive (sterile or not) 
neutrino. Therefore, such a test is particularly powerful 
and worth addressing as the next step of our analysis. 

It is worth noting that help investigating the viability 
of our MONDian vector models may come from a non 
astrophysical experiment. Indeed, one of the key ingre- 
dients in both MONDian vector models is the presence 
of massive neutrinos with m v ~ 2 eV. The KATRIN 
experiment [57j on the tritium (3 decay should be able 



to constrain the electron neutrino mass with a sensitity 
of ~ 0.2 eV. Should this experiment indeed find that 
neutrinos are less massive than, e.g., ~ 1 eV, our model 
could be in serious trouble unless one assumes that a 
single massive sterile neutrino does indeed exist. 

As a final remark, one could note that we have titled 
our paper with a question so that the reader could now 
claim for an answer. Unfortunately, because of the un- 
certainties on the use of the distance priors and the small 
redshift range probed, the unique answer we can give af- 
ter our analysis is only a (somewhat frustrating) maybe. 
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